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Congestion in systems with one server, exemplified by simple queues, in- 
dividual telephone trunks and particle counters, is considered without any 
restrictions on the statistical character of the incoming load. Elementary 
methods establish formulas and equations describing probabilities of delay and 
loss. These methods deemphasize special statistical models and yield a gen- 
eral theory. In spite of this generality, it is attempted to give intuitive proofs 
and extensive explanations of the physical significance of formulas, as well 
as rigorous derivations. The theory is applied to specific models to obtain 
illustrative new residts. 

I. INTRODUCTION AND SUMMARY 

Congestion theory is the study of mathematical models of service 
systems such as telephone central offices, waiting lines and trunk groups. 
It has two practical uses: to provide engineers with specific mathemati- 
cal results — curves and tables — for designing actual systems, and to 
provide a general framework of concepts into which new problems can 
be fitted and in which current problems can be solved. Corresponding 
to these two uses are two kinds of results: specific results pertaining to 
special models and general theorems valid for many models. 

Most of the present literature of congestion theory consists of specific 
results resting on particular statistical assumptions about the traffic in 
the service system under study. Indeed, few results are known in con- 
gestion theory that do not depend on special statistical assumptions, 
such as negative exponential distributions or independent random vari- 
ables. In this paper we obtain some mathematical results that are free 
of statistical restrictions, but which apply only to a limited class of sys- 
tems. These results concern general stochastic processes in traffic sys- 
tems, such as counters and simple queues, that have only one server. 

The limitation to a single server is severe, and may be difficult to 
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overcome. We hope that appropriate analogs of the methods used here, 
applied to stochastic processes in several dimensions, will yield statisti- 
cally general theories for systems with many servers. Some systems, for 
example, queues with many servers in parallel, cannot be given a theory 
free of special assumptions without the use of multidimensional proc- 
esses; other systems have a structure that permits us, in principle, to 
treat them by "iterating" results for one server. In Section VIII we 
suggest, as a conjecture, how the general theory we give for one tele- 
phone trunk might be applied several times to give a general loss formula 
for a group of trunks. 

The classical problems of telephone traffic engineering usually involve 
many servers, and so theories dealing with single servers have been of 
peripheral interest to traffic engineering in the past. However, modern 
increases in the speed of components have made possible electronic tele- 
phone exchanges, now under study and development. By dint of their 
speed, some of these systems either have a single active server (the 
"common control"), or may reasonably be viewed as a set of (relatively 
independent) single servers. Therefore, the case of one server is in fact 
of immediate and practical interest. 

The aims of this paper are three: (a) to describe a new general ap- 
proach to certain congestion problems; (b) to show that this approach, 
although general, can nevertheless be presented in a relatively elemen- 
tary way that makes it available to many persons; (c) to illustrate how 
the new approach yields specific results, both new and known. What 
follows is written only partly as a contribution to the mathematical 
analysis of congestion. It is also a frankly tutorial account, aimed at 
increasing the public understanding of congestion by first steering atten- 
tion away from special statistical models and then obtaining a general 
theory. Such a point of view, it is hoped, will yield new methods in prob- 
lems other than congestion. 

When a general theory can be given, it should be useful in several 
ways: in increasing our understanding of complex systems; in obtaining 
new specific results, curves, tables, etc; in extending theory to cover 
interesting cases that are known to be inadequately described by exist- 
ing results. At first acquaintance, the theorems of such a general theory 
may not resemble "results" at all; that is, they may not seem to be 
facts that one could obviously and easily use to solve a real problem. A 
general theory is really a tool or principle expressing the essence or 
structure of a system; properly explained and used, this tool will yield 
formulas and other specifics with which problems can be treated. 

Let us give examples of questions we shall be able to answer and re- 
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suits we shall obtain. As a first example of increased understanding, we 
shall be able to answer, affirmatively, the question 

Is there a small number of specific features of the incoming 
traffic which suffices, quite generally, to determine delay dis- 
tributions and loss probabilities? 

We shall exhibit these features, and give delay and loss probabilities in 

terms of them. 

As a second example, of both increased understanding and a new case, 
let us consider the matter of Markov stochastic processes. A Markov 
process has a particularly simple probabilistic structure, in that all in- 
formation about the past history of the randomly moving point is irrele- 
vant to its future development, if its present position is known. Such 
processes have been very useful in congestion theory, because of the 
convenient functional equations (such as "statistical equilibrium" equa- 
tions, continuity equations, renewal equations) that are associated with 
them. However, it is also true that congestion theory has occupied itself 
almost exclusively with Markov processes, or with the generalizations of 
these called regenerative, and semi-Markov, processes. Thus the follow- 
ing question has doubtless been asked: 

What functional equations can be derived for non-Markov 
problems in delay, or in loss? 

In answer, we show that the probability that a queue be empty at t 
satisfies, quite generally, a Volterra equation of the first kind. This 
equation coincides with a known equation in the Markov case, as it 
should; similar results hold for loss. 

By way of specific results, we prove that, in many cases of practical 
interest described by Markov processes, the solution of the Volterra 
equation has the intuitive form 

Pr{ queue is empty at t) = 

i L total load in [0,0 \ 
average of the greater of and (1 I 

+ initial load (chance that total load in [0,0 is at most t). 

In these cases, it is no longer necessary to solve the Volterra equation; 
one merely computes the quantities on the right-hand side of the for- 
mula just given. The "initial load" is the time the queue would take to 
empty if no customers arrived after time zero; the "total load in [0,0" 
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is the initial load plus the sum of all the service-times of customers 
arriving in (0,/). 

This result is illustrated by Fig. 1, which shows the probability that 
a queue is empty as a function of time for two well-known cases, nega- 
tive exponential service times and constant service times. In each case, 
arrivals are in a Poisson process at the rate X, and the mean service 
time is b. The curves suggest that the approach to equilibrium is con- 
siderably faster with constant service times than with exponential serv- 
ice times. The "exponential" curve was computed from (48), and the 
"constant" curve from the formula given above. Both curves approach 
the same limiting value, 0.750. 

The physical background of intended applications is discussed briefly 
in Section II, while Section III describes the cumulative load or traffic. 
Sections IV and V are devoted, respectively, to equations characterizing 
delay and loss operation in terms of the offered load; these sections are 
purely descriptive, and no probability is involved. Probability first enters 
in Section VI, where we discuss the general nature of the problems we 
try to solve, the methods we use and the results we obtain. Formulas 
and equations are stated and explained for probabilities of delay in 
Section VII, and for probabilities arising in loss operation in Section 
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Fig. 1 — Probability that a queue, starting empty at 0, is empty at t, as a func- 
tion of t for negative exponential and constant service times, with A = 1, b = 0.25. 
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VIII. Proofs of, and additional comments on, the results for delay appear 
in Section IX, and a number of specific formulas are deduced in Section 
X. Precise derivations of the results for loss are in Section XI, with an 
additional limit theorem being presented in Section XII. 

II. PHYSICAL BACKGROUND 

We shall discuss results applicable to the following type of system: 
there is given a machine that is either idle or busy; every so often, at 
random instants, someone tries to use the machine; and the machine can 
be operated in two ways, called loss and delay. In delay operation, people 
who find the machine busy wait for their turns in order of arrival, and 
then make it busy (by using it) for the originally intended length of 
time; in loss operation, people who find the machine busy are sent away 
or "lost," while those who find the machine idle seize it and use it for 
a random time, the service time. 

The system to be studied has been described in terms of people who 
try to use a machine, but the same structure appears in other applica- 
tions. In telephony, for example, the machine to be used may be a 
marker, and the "people" arriving may be registers filled with dialed 
digits; telephone engineers usually refer to the service time as the hold- 
ing time. In the study of radioactivity, the machine is an ionization 
counter, the "people" are impinging particles, and the service time 
(called the dead time of the counter) is a length of time after the register- 
ing of a particle during which the counter cannot count, so that it misses 
particles arriving during this interval. 

The quantities of engineering interest in these models vary with the 
application, but the following ones seem to be of general importance: for 
loss operation, the chance of loss, i.e., the probability that an arriving 
customer find the machine busy, and the probability distribution of the 
busy period, i.e., the amount of time that must elapse until the machine 
next becomes idle; for delay operation, the distribution of waiting time. 

III. THE CUMULATIVE LOAD, K(t) 

Before we can study the two modes of operation, loss and delay, we 
must describe the offered load of work, or the arriving traffic. This is 
done most easily by using a step function, K(t), which is nondecreasing 
and left-continuous. The locations of the jumps are the epochs of arrival 
of customers, and the magnitudes of the jumps are the service times or 
the lengths of time the machine is made busy. Equivalently, the offered 
load is completely determined by the arrival epoch t k and the service- 
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time Sk of the kth. arriving unit, k = 1,2, • • • . This situation is depicted 
in Fig. 2. We shall use K(t) to describe the load for both loss and delay 
operation. 



IV. INTEGRAL EQUATION DESCRIBING DELAY OPERATION 

As a mathematical description of the delays to be encountered under 
delay operation, we use the virtual waiting time, Wit), which can be 
defined as the time a unit would have to wait for "service" if it arrived 
at time t. At the epoch t„ of arrival of the nth unit, Wit) jumps upward 
discontinuously an amount equal to S„ , the work or service time of the 
unit. Otherwise, W{t) has slope —1 if it is positive; if it reaches zero, 
it stays equal to zero until the next jump of the load function Kit). 
Corresponding graphs of Kit) and W(t) appear in Fig. 3. 

If Kit) is interpreted as the work offered in the interval (0,0 , then 
Wit) can be thought of as the amount of work remaining to be done 
at time t. In terms of this interpretation, it can be seen that 

work remaining at t = total work load offered up to t — elapsed time 

+ total time during which machine was idle in i0,t). 

The machine is idle when, and only when, Win) — 0. Then, formally, 
Wit) is defined in terms of Kit) by the integral equation 



W(t) = Kit) -t+ I U[-Wiu)]du, t^O, 
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Fig. 2 — Cumulative load. 
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Fig. 3 — (a) Graph of K(t) vs. t; (b) graph of W(t) vs. t. 

where U(t) is the unit step function; i.e., U(x) = 1 for x ^ 0, and 
U(x) = otherwise. For simplicity, we have set TP(0) = K(0). 

It has been shown by Reich 1 that the solution of (1) is, if K(x) — x 
has a zero in (0,0, 

W(t) = max [K(t) - K(x) - t + x]. 

This fact may be interpreted physically as follows: the quantity in 
brackets, [K(t) — K(x) — t + x], is, if positive, the excess of arriving 
load in the interval (x,t) over the elapsed time t — x;\t is therefore the 
overload in (x,t). Reich's formula then says essentially that 

delay at t = max [overload in (x,t)]. 

The relationship between the waiting time W(t) and the offered 
traffic K(t) can be further elucidated graphically by reference to Fig. 4. 
The light solid line shows K(t) — t, the traffic offered up to time t minus 
the traffic that could have been served if the server had been kept busy 
throughout the interval (0,0- It is supposed in Fig. 4 that the server 
starts busy at t = 0. It is busy until t = a. At this point, the server be- 
comes idle and K(t) — t turns negative, and its negative value is the 
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Fig. 4 — Relationship between waiting time and offered traffic. 



negative of the idle time. At t = b, more traffic is offered and K(t) — t 
jumps up. 

The heavy dashed line represents the waiting time at t, W(t); W(t) 
can also be thought of as the traffic in hand and yet to be served. It 
can never go negative. It is equal to K(t) — t before a, and is zero from 
a to b. At b, it jumps up, remaining above and parallel to K(t) — t 
until t = d, when the server goes idle again. At b, W(t) is above K(t) — 
t by exactly the amount by which K(t) — t was most negative. At d, 
when W(t) reaches zero, K(t) — t is just reaching its previous minimum, 
and K(d) - d = K(b) - b. 

During the interval (d,e), W(t) remains at zero, while K(t) — t be- 
comes more negative, establishing new minima as it goes and building 
up more idle time. At t = e, K(t) — t and W(t) both jump up; W(t) 
is again parallel to K{t) — t, but it is now above it by an amount equal 
to the negative of the last minimum, K(e) — e. 

In Fig. 4, 

min [K(x) — x], t ^ a, 

is shown as a light dashed line. It is a monotone, nonincreasing function 
of t, and is the negative of the total idle time up to time t. To account 
for the period t < a, when K(t) — t has not yet become negative and 
the server has not yet been idle, we write 

W(t) = K(t) - t - min {0, min [K(x) - x]}, 

0<x<t 
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and thus obtain another representation for the delay, i.e., a solution of 

(1). 
In a manner similar to that of Fig. 4, Fig. 5 depicts, simultaneously, 

the offered load K(t) in a light solid line; the waiting-time W(t) in a 
heavy solid line; the negative of the accumulated idle time in a heavy 
dashed line; and the "load-time excess" K(t) - t, when it does not 
coincide with the negative of the idle time, in a light dashed line. The 
terminology in Fig. 5 has been purposely chosen to suggest an interpre- 
tation in terms of inventory or storage theory: W(t) is the real backlog 
(of orders, say,), K(t) is the cumulative amount ordered and K(i) — t 
might be termed the load time excess or the virtual backlog. Then 

real backlog = virtual backlog + accumulated idle time. 

V. INTEGRAL EQUATION DESCRIBING LOSS OPERATION 

For loss operation, we use A(t), the service time remaining at t, as 
an indicator of the condition of the machine. We may define A(t) as 
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Fig. 5 — Relationships among offered load, waiting time, negative of accumu- 
lated idle time and "load time excess". 
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the amount of time remaining until the machine next becomes idle, 
with the convention that A(t) = means that the machine is idle at 
t. At epochs of arrivals, A(t) jumps upward discontinuously if the arrival 
in question finds the machine idle, the amount of jump being the service 
time of the arriving customer. Arrivals that find the machine busy do 
not affect A(t). While A(t) is positive, it has slope — 1; when it reaches 
zero, it stays equal to zero until the next jump. A graph of A (t) consists 
of isosceles "right triangles" laid on their sides; corresponding graphs 
of K(t) and A(t) appear in Fig. 6. For the study of A(t), it is more 
convenient to define K(t) to be right-continuous. 

In terms of the interpretation of K(t) as cumulative work load it is 
easy to see that 

work remaining to be done at t = total work load offered up to t 

— load missed in (0,0 because machine was busy 

— elapsed time + total time during which machine was idle in (0,0- 
This means that A(t) is defined in terms of K(t) by the equation 

A(t) = K(t) - [ {1 - U[-A(u)]} dK(u) - t 

(2) 



+ [ U[-A(u)]du, 

Jo 



where, as before, U(t) is the unit step function, and we have set ^4(0) = 
7i(0). 

VI. CHARACTER OF THE GENERAL RESULTS 

Models of waiting lines and telephone traffic usually contain explicit 
assumptions about the statistical nature of the offered load K(t). For 
instance, Erlang's original models 2 amount to assuming that the inter- 
arrival times (t„ — t n -i) are all independent with the same negative 
exponential distribution; and similarly for the service times. These as- 
sumptions give a class of models whose parameters are the means of 
the negative exponential distributions. 

A broader class of models is specified by retaining the assumptions 
that the interarrival times be independent and identically distributed 
(and similarly for service times), but allowing any distribution, not 
just the negative exponential. The interarrival and service time distri- 
butions may still be said to "parametrize" this broader class of models, 
since their choice determines a model in the class. 
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Fig. 6 — Graphs of (a) K(t) vs. t and (b) W(t) vs. t for loss operation. 

As more general assumptions even than those of the last paragraph 
are considered, it becomes extremely laborious first to specify the model 
and then to compute interesting quantities such as distributions of delay, 
probabilities of loss, etc. So, instead of looking for ways of exactly char- 
acterizing the model, we can try to search directly for simple ways of 
expressing the quantities of interest in terms of the model. Since the 
probabilities 



¥i-\W(t) ^ w) 
Pr{A(t) ^ w] 



(3) 
(4) 



are what we actually wish to compute from the model, the question 
arises whether this calculation can be made without first specifying the 
entire probabilistic structure of K(t). The following intuitive argument 
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can be adduced for answering "yes". Both W(t) and A(t) are denned 
in terms of the load by very special relationships that are expressed in 
the integral equations we have given. Hence, no matter what are the 
statistical features of K(t), it is likely that the distributions of W(t) and 
A(i) depend only on some very particular, physically interpretable 
statistical functions associated with K(t). It is not obvious that such 
an economy can be made in the generality we desire. 

Two principal results, described in the next sections, state that the 
probabilities (3) and (4) can in fact be given fairly simple expressions 
that are generally valid for any load process. These expressions depend 
only on certain special functions obtainable from the statistical structure 
of the load K(t). Each function has a definite intuitive or physical signifi- 
cance, to be given later. These statistical functions achieve the desired 
economy of description, because we can state that the desired probabili- 
ties depend only on the features oi K(t) expressed in the functions. For 
the purposes of calculating (3) and (4), we do not need the entire proba- 
bilistic structure of K(t), but only the small relevant part. 

Another way of putting the problem that we have attempted to 
solve is: for general load processes K(t), what small amount of informa- 
tion about the statistical nature of K(t) will determine ~Pr{W(t) ^ w) 
for all t and w, and how is this computation to be made? [A similar 
question arises for A(t).] The statistical functions are then the informa- 
tion about K(t); the formulas for the probabilities (3) and (4) indicate 
the method of computation. Specifically, (1) defines a transformation of 
a stochastic process, K(t), of service times and arrival epochs into an- 
other stochastic process, W(t), of waiting times. For each t, there is an 
operator or formula that gives the distribution of W(t) in terms of suit- 
able fundamental statistical functions associated with K(u) for u S t. 
The principal problem is to find the form of the operator and the char- 
acter of the fundamental functions. The answer to this problem should 
depend only on the integral equation (1) and on the fact that K(t) is a 
nondecreasing step function. It should depend on no special features of 
the probability measure for K(t) except those implied by this last prop- 
erty. Accordingly, we shall assume at first that K(t) is a random, nonde- 
creasing step function; its only statistical peculiarity is that it is a non- 
decreasing step function. 

There already exists an extensive literature dealing with probabilistic 
models of telephone trunks, electronic counters and similar machines; 
we shall attempt to relate our approach to this literature. The models 
used to date have usually included strong hypotheses of independence 
or else used special distribution functions, such as the negative exponen- 
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tial. The results obtained depend methodologically, almost without ex- 
ception, on the possibility of writing Kolmogorov equations for a Markov 
process, or integral equations of renewal type for probabilities or expec- 
tations associated with a regenerative process (Smith ). The extensive 
work of Tak&cs 4 ' 6 ' 6 falls entirely in this category (see also the references 
therein). 

Our approach is, in a sense, an inversion of the usual method described 
above. The latter consists in first doing probability theory to set up 
Kolmogorov or renewal equations, and then doing analysis to solve the 
equations. We can, however, achieve greater generality by taking maxi- 
mum advantage of the fact that the processes of interest, W(t) and A(t), 
already satisfy (1) and (2), respectively. This we do by careful analysis 
of W(t) and ^4(0 themselves, by in effect performing some of our analy- 
sis in the domain of random variables, and taking averages only at 
convenient points. 

VII. PROBABILITIES OF DELAY 

It has been shown 7 that Pr{W(t) ^ 10} can be expressed in terms of 
the two statistical functions 

Fr{K(t) ^ w), (5) 

R(t,u,w) = Vx\K(t) - K(u) - t + u S w | W(u) = 0|. ((>) 

We first clarify the physical significance of the functions (5) and (G) 
entering into the representation of Pr{TF(/) ^ w}. The first is simply 
the probability distribution of the cumulative load, K(t). The second, 
written as 

R{t,u,w) = Vv{[K(t) - (} - [K(u) - u] ^ w\W(u) - 0} 

is the (conditional) probability distribution of the change in the "over- 
load" that occurs in the interval (u,t), given that the server was idle at 
time u. Here we have used the term "overload" to mean the amount of 
work that would have been left still undone even if the traffic had been 
so spaced that there was no idle time. This overload is represented by 
K(t) — t. A negative value of the overload represents the negative of the 
idle time that would have occurred if the traffic spacing had permitted 
the server to finish all of the work by time t. 

The significance of R(t,u,w) can be made more apparent by describing 
how to measure it. We look at many identical copies of the system at a 
time u when the server is idle, and we measure the totaled service times 
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of all customers arriving in the next interval (u,t), fort>u; the fraction 
of these numbers that is in the range (0, t — u + w) is R(t,u,w). 
The delay W(t) is never negative, and so 

Vv{W(t) ^ w) = 0, if w < 0. (7) 

For w > 0, Pr{TT(0 ^ u>} is given by 

Pr{TT(0 ^w) = Pr{K(0 - * ^ w) 

(8) 



- — f R(t,u,w) Pr{W(u) = 0} dw. 

aw H 



The chance, Pr{JF(w) = 0}, that the server will be idle at time u satis- 
fies the Volterra equation of the first kind, for — t 2£ w 2g 0, 

average of max [0,w — X(<) + £] = 

.t+w (9) 

R(t,u,w,) Pt{W(u) = 0} dw, 



/ 



and the left side of (9) is expressible in terms of (5) as 



i+i. 



#{max [0,w - K(t) + i\) = \ Fr[K(t) ^ u} du, 

Jo 

where E\u\ is the expectation of u. 

Once the "basic" functions (5) and (6) are known, the computation 
of Pr{W(t) ^ w\ proceeds by first solving the integral equation (9) 
for the chance Ft{W(u) = 0} that the server will be idle at time u. This 
probability can then be used in (8) to give Fr{W(t) ^ w). 

Proofs for, and extensive explanations of, these results for delay have 
been deferred until Section IX, while we continue by discussing results 
for loss operation. The proofs to be given are new, and much simpler 
than those of Ref. 7. They also make it easier to exhibit the physical 
significance of the formulas and functions arising. 

VIII. FORMULAS FOR LOSS OPERATION 

We shall show that Pr{A(t) ^ w) can be expressed in terms of two 
kernels, R(t,u) and Q(t,u), as follows: 

Pr{A(«) ^ w) = Pr{A(0) ^ t + w] 

r* (10) 

- [1 -R(t+ w,u)]dE{S(u)\, 

Jo 

where E{S(u)} is the average number of units that arrive in (0,w] and 
find the machine idle. It satisfies the integral equation 
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E[8(t)} = Prfo ^ t\ + [ Q(t,u)dE{S(u)}, (11) 

where 2/1 is the epoch of the first successful arrival. 

The kernel R(t,u) may be interpreted as a rigorous version of 

Pr{ service time of a successful unit arriving at u is ^ / — u | a 
successful unit has arrived at u] , 

and Q(t,u) may be thought of as a rigorous version of 

Pr{ next successful arrival after u occurs before t | a success- 
ful unit has arrived at u) . 

The first term on the right in (10) is self-explanatory. Precise definitions 
of R( • , •) and Q( ■ , ■ ) are given in Section XI. 

To explain (10) itself, we observe that A(0) > t + w implies .4(0 = 
.4(0) — t > w, and hence 

Fr{A(t) ^w}^ Pr{A(0) ^ t + w). 

The integral term in (10) is therefore a correction to the overestimate, 
Pr{i4(0) ^ t + w). From the interpretation of R(t,u), and that of 
dE{S(u)\ as the "density" of successful arrivals at u, we see that (10) 
can be rendered in words as 

Prjwork time left at t ^ w\ = Prfwork time left at ^ t + w} 

-Prfsomc successful arrival during (0,t] stays beyond time t + w). 

It is reasonable to suspect that, if the load process K(t ) has some weak 
stationarity properties, and if certain "averages" exist, then 

Pi{A(t) ^ w) 

has, for each w ^ 0, a nonzero limit as t — > 00 . For w = 0, one expects 
intuitively that this limit will have the form 

Pr{A(a.) = 0} = 



+ 6' 



where a,b are constants which can be interpreted as follows : if we watch 
the process A(t), we notice that periods of time during which A(t) = 
alternate with those during which A(t) > 0; then a is the average 
length of a period during which A(t) = 0, and b is the average length 
of a period during which .4(0 exceeds 0. 

These conjectures are justified in Section XII, where it is shown that, 
if the kernels R(t,u) and Q(t,u) used in (10) and (11) are functions of 
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(t — u) only, and if 

1 - R(x) (12) 

is integrable over (0, oo ) , then 

limPr{A(0 £ w} 

i-*oo 

exists and has the form 

a + I [1 - R(u)] du 
Jq 

a + b 
with 



b = f [1 - R(u)]du, 

Jo 

a+b = / [1 - Q(u)]du. 
Jo 



We have already pointed out that R(t,u) can be interpreted as 

Pr{ service time of a successful unit which arrives at u is (13) 
^ t — u | a successful unit arrived at u} . 

Then R(t,u) = R(l — u) states that (13) does not depend on the first 
and third occurrences of u therein; i.e., that service times have the same 
distribution no matter when they begin. Thus, the dependence of R(t,u) 
on (t — u) only is a weak sort of stationarity property, and R(x) can 
be interpreted as the probability distribution of service times of success- 
ful units. Then 

6 = / [1 — R{u)] du = {average service time of successful units} , 
Jo 

since when the mean of a positive variate exists it equals the integral of 
the "tail" of its distribution. 

In a similar way, the fact that Q(t,u) is a difference kernel may be 
interpreted as a stationarity condition, and Q(x) can be thought of as 
the distribution function of the intervals between successful arrivals, so 
that 

a+b = [ [1 - Q(x)] dx 
Jo 

= {average interval between successful arrivals} 
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and 

a = average length of an idle period. 
The preceding discussion suggests that we use 

chance of loss = — ; — ; (14) 

a + b 

as a natural measure of the probability of loss for the single telephone 
trunk or particle counter that we are considering. Of course, any engi- 
neer would have used (14) to describe loss, justifying it by intuitive 
arguments. This fact does not detract from our result, which gives some 
idea of the weak assumptions that are sufficient for proving (14). 

Suppose that, instead of having only one machine, we had N ^ 1 
machines, and used them in a fixed serial order of preference. That is, 
arrivals finding the first n machines busy try the (n + l)th. Such a 
situation arises in telephony, for instance: the "machines" are telephone 
trunk lines, the "arriving units" are attempts to place a call and the 
work or service times are the holding times of calls. Each trunk is then 
receiving the traffic overflowing the previous trunks in the ordering. 
Our theory then applies to each trunk considered by itself, and if the 
conditions for the validity of (14) obtain for each trunk, the chance of 
loss for the whole group must have the form 

n ' 



kha„ + &»' 

where 

__ f average service time of units that find the first (n — 1) 
n ~ \trunks busy, the nth idle 

_ /average time interval between calls accommodated on\ 
a„ + b n - \ thenthtrunk J' 

Formulas (10) and (11) have been essentially proved 8 by involved 
arguments using the integral equation (2) defining A(t). We shall give 
a simple heuristic derivation in this section and a rigorous one in Sec- 
tion XI. 

In order to explain the results, we first recall that the process A{t) 
consists of alternating intervals during which A(t) is first zero, then 
positive with slope —1, then zero again, and so on. The next arrival to 
find the machine idle always makes it busy again. 

We are interested in expressing Pr{A(t) ^ w], and so we search for 
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other ways of specifying the event \A(t) ^ w). We first consider those 
cases in which ^1(0) = 0; that is, the system starts empty, as in Fig. 7. 
Then it is not hard to see that A(t) will be less than or equal to w only 
if the number S(t) of successful arrivals during (0,2] equals the number 
of successful arrivals in (0,t] that have left the system by time t + w. 
In other words, if the machine is idle at t = 0, then the work A(t) re- 
maining at time t is less than or equal to w if and only if all people who 
arrived to find it idle in (0,i\ are finished with it by time t + w. We shall 
set 

d(t,t + w) = number of successful arrivals in (0,1] who have 
left the system by time t + w. 



Then, if 4(0) = 0, the events 
[A(t) ^ w\ and 



S(t) = 6(t,t + w)} 



(15) 



are the same. 

If A(0) > 0, the system starts busy, and the graph of A(t) appears 
as in Fig. 8. Assume first that A(0) > t; then, of course, A(t) — 
— t + ^4.(0), because the machine has been busy since t = 0, and is not 
yet finished. If t ^ A(0), though, the machine became idle at t = A(0). 
In the first instance, S(t) = d(t,t + w) = 0, because there have not 
yet been any successful arrivals, and A(t) ^ w if and only if 

t < A(0) ^ t + w. 

In the second instance, the argument we used for the case A(0) = 
applies, and (15) holds. We now average these cases according to their 




y. * 2 y 2 

TIME, t — *• 

Fig. 7 — Graph of A(t) vs. t when system starts empty. 
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TIME, t — *• 

Fig. 8 — Graph of .4(0 vs. I when system starts busy. 

probabilities; since either S(t) — 6(t,t + w) = or S(t) — 6(t, 
t + w) = 1, we find 

Vr{A(t) ^w) =Vr{t<A(0) ^t + w and S(t) = 6(t,t + w)\ ' 

(16) 
4- Pr{A(0) ^ f and S(t) = d(t,t + w)}. 

Where { ■ • •} is an event, let x! • • •} be its characteristic function, de- 
fined as 1 if the event happens, and otherwise. Then the first term on 
the right in (16) is 

average of (x{t < A(0) ^ t + w} -[1 - S(t) + d(t,t 4- w)]), 

and the second is 

average of (x{A(0) ^ t] -[1 - S(t) + d(t,t + w)]), 

because 

1 - S(t) 4- e(l,t + w) = if 5(0 ^ 0(M 4- w). 

Therefore 

Pr{X(«) ^ w) = Pr{A(0) ^ t + w\ - #{«(*)} + W(M + «01- (17) 

Thus, we have expressed Pr{ A(t) ^ w) in terms of the initial distribu- 
tion and the average values of S(t) and 6(t,t + w). 

To prove (10), it remains to express the average of 6(t,t -f- w) as the 
integral of a kernel with respect to the average of S(t). To do this, let 
u be a point in (0,t\, and suppose that a successful arrival occurs at u. 
Such an arrival can only contribute to the average of 6(t,t + w) if it 
leaves before / + w. The proportion of such arrivals is just 
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Pr{ service time of a successful unit which arrives at u is 

^ t + w — u I a successful arrival occurred at u) (18) 

= R(t + w,u). 

The "density" of successful arrivals at u is dE[S(u)} ; therefore, 

#{0(M + w)) = [ R(t + w,u) dE{S(u)}, 
Jo 

which proves (10). 

We continue with a heuristic derivation of the integral equation (11) 
for E[S(t)}. First we notice that the event {A(t) = 0} can occur in two 
ways: either some successful arrival has occurred in (0,t], or none has. 
If none has, then either A(-) started idle at and is still idle at t, or 
else it started busy at 0, became idle at the point .4(0) < t and is still 
idle at t. If h is the first arrival in (0, ■» ) and yi the first successful arrival, 
the chance that no successful arrivals occurred in (0,t] and A(t) = is 

Pr{4(0) = and t % > t) + Pr{0 < A(0) £ t and y x > t). 

Assuming some successful arrival did occur in (0,/], suppose it occurred 
at u. Such an arrival is only relevant to the event {A(t) = 0} if it is the 
last such arrival in (0,i\, and if the service time of the customer then 
arriving is at most (t — u). The proportion of such "relevant" arrivals 

is 

Pr{ service time of successful arrival occurring at u is ^ t — u 
and no more customers arrive in the time interval between , _. 
his departure and t \ a successful arrival occurred at u) 
= G(t,u). 

As before, we now argue that the density of successful arrivals at u is 
dE{S(u)} and so, using (19) and noting that yi = h if A(0) = 0, we 
find 

Pr{A(0 = 0} = Vr{A(0) ^ t, Vl > t) + [ G(t,u) dE{S(u)}. (20) 

Jo 

By combining this result with (10) for w = 0, we obtain an integral 
equation for E{S(u)} : 

E{S(t)} = Pr{2/i ^ t} + [ Q(t,u) dE{S(u)}, (21) 

Jo 

where the kernel Q(t,u) = R(t,u) — G(t,u). By examining the inter- 
pretations (18) and (20) of R(t,u) and G(t,u), it can be seen that 
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Q(t,u) = Prjncxt successful arrival after u occurs before / | a 
successful arrival occurred at u) . 

Hence, Q(t,u) should be a distribution function in t. Proof of this, to- 
gether with discussions of (18) and (20), appears in Section XI. 

IX. PROOF AND DISCUSSION OF THE RESULTS FOR DELAY 

The proof and explanation of (8) and (9) depend on two simple pre- 
liminary results. The first of these is as follows: let x be any nonnegative 
random variable; then, for ij ^ 0, 

t Pr(.r g u) du = #{max (0,y - x)}. (22) 

Jo 

This formula states that the area under the (cumulative) distribution of 
x to the left of y is just the average value of the greater of zero and 
y — x. This is easily seen from an integration by parts: 

ru \y rv 

I Pr{z ^ u\ du = u Fr{x ^ u}\ — I ud Pr{as ^ u} . 

Jo |o «<0 

To begin the proof we note that the total idle time T(t) represented by 
the term 



T{t) = f U[-W(u)]du 
Jo 



in the integral equation ( 1 ) , is always nonnegative, so that 

Wit) ^ K(t) - t, (23) 

and also 

Pr{TT'(0 ^ w] ^ Fr{K(t) - t ^ w}. 

Now in Fig. 9 the larger area represents the event { K(t) — t ^ w) and 
the smaller one inside it represents the event {W(<) ^ w). This latter 
event is included in the former because, if W(t) ^ w, then 

K(t) - t ^ iv, 

by (23). The difference between the two areas represents the event 

\K(t) - t ^ w < K(t) - t + T(t)}, 

and so 

Fr[W(t) ^w}= Pr{K(t) - t ^ w) 

(24) 
- Pr \K(t) - t ^ w < K(t) - t + T(t)}. 



(25) 
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We next observe that (8) is the derivative with respect to w of 
/ Pr[W(t) ^ u) du = \ Pr{K(t) - t ^ u] du 

Jo Jo 

- { R(t,u,w) Pr{W(u) = 0} du, 
Jo 

which may be written, using (24) and taking the condition inside, as 

E{max[0,w - K(t) + t - T(t)] - max[0,u; - K(t) + l]\ 

= - f Pr{A'(/) - K(u) - t + u ^ w and W(u) = 0} du. (26 ' ) 

Jo 

Thus (8) and (9) are established if we can prove (26). 

To establish (26) we need the second preliminary result, a general 
property of monotone continuous functions. 

Lemma: If F(t) is continuous and monotone increasing and F(0) = 0, 
then, for any x ^ and t ^ 0, 

max[0,.r - F(t)) = x - I U[x - F(y)] dF(y), (27) 

Jo 

where U(y) = 1 for y ^ 0, and U(y) = for y < 0. 
Proof: We note that, as y increases, the integrand in (27) is unity until 
either x = F(y) or y = I, whichever occurs first, and it is zero thereafter. 
If x = F(y) occurs first, then the integral equals x; if y = t occurs first, 




Fig. 9 — Graphical representation of events \K(t) -Jg»] and j W(t) — t ^ w\. 
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the integral equals F(t). Hence, 

( U[x - F(y)] dF(y) = mm[x,F(t)]. (28) 

Jo 

For x ^ it can be seen that 

mm[x,F(t)] = x - max[0,.r - Fit)], 

and this proves the lemma. 

To use this result in proving (2G) we interpret F(t) in the lemma as 
T(t), the total idle time during (0,0, which is a continuous increasing 
function, with 7'(0) = 0. We next consider the expressions, for w ^ — t, 

A = w - Kit) + t - [ U[w - Kit) + t - T(u)] dT(u), 
Jo 

B = w - Kit) + t + T(t) 

- { U[w - K(t) + / + r(0 - Tiu)]dTiu). 
Jo 

If w — Kit) + / ^ 0, then, by the lemma, 

A - B = max[0,w - Kit) + t - T(t)) - max[0,w - Kit) + t]. (29) 

However, since w — Kit) + t ^ 0, we see that the integrand in B is 
always unity, so that B equals w — K (t) + t, and 

A - B - -f U[w - Kit) + t - Tiu)] dTiu). 

Jo 
Hence w — Kit) + t ^ implies 
max[0,w - Kit) + t - Tit)] - ma\[0,w - Kit) + t] 

(30) 



= - f U[w - Kit) + t - Tiu)] dTiu). 
Jq 



However, if w — Kit) + /• < 0, then both sides of (30) vanish and so 
(30) holds generally, for any w — Kit) + t. 
From (1) it is evident that, while W{u) is zero, then 

T{u) = u - Kiu); 

so the integral in (30) is 

f U[w - Kit) + t + Kiu) -u]^p^du, (31) 

Jo du 
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whore dT(u)/du is one if W(u) = 0, and is zero if not; i.e., 

( Ej^l = U[-W(u)}. (32) 

du 

We recall that, if x and y are random variables, then 

Pr{.r ^ X, and y ^ X 2 | = E{U[h - x]U[M - y}}; (33) 

i.e., the joint distribution is the average of the product 

C7[X, - x]-U[X,- y]. 

If we now average (30) and bear in mind (31), (32) and (33), we ob- 
tain (26) for all w ^ — t. For negative w, 

E[max[0,w - K(t) + t - T(t)]\ = 0, 

so that formula (26) takes the form, for — t ^ w ^ 0, 

#{max[0,w - K(t) + t]} = 

rt+ W (34) 

Pr{A(0 - K{u) - t + u ^w and W(u) = 0} du. 
Jo 

This is (9), and we have completed the proof of the results stated at the 
beginning of this section. 

It can be seen from (8) and (24), and from Fig. 9, that 

— f Pr(A(0 - K(u) - t + u ^w and W(u) = 0} du (35) 

dW Jo 

is the correction term to the overestimate Pr{K(t) — t ^ w) forPr{TT(<) 
^ w). [See (24).] It is the probability of the event 

\K(t) -t£w < K(t) - t + T(t)}, 

represented by the difference between the areas in Fig. 9. Now, the 
presence of the derivative d/dw in (24) is explained by the fact that 
(24) is the derivative of (25) for w ^ 0, and thus is due to our use of 
(22). However, it is not obvious intuitively why, in (35), the rest of the 
term after the d/dw, should be a time integral. 

An explanation of this can be obtained from (26), which expresses 
the average of the random variable 

a = max[0,u> - K{t) + t - T{t)\ - max[0,w - K(t) + *]. 

It can be seen that 
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-7X0 
a = <j - w + #(<) - i 





if w> K(«) - t + T(0 = TT(0 
if K(t) - t <w < Wit) 
if w < K(t) - t. 



A graph of a as a function of w is shown in Fig. 10. The point K(t) — t 
at which a starts downward may, of course, be negative, although it is 
positive in the figure. 

Thus, a is a negative quantity whose magnitude is no greater than 
T(t), the total idle time prior to t. Now, if a were in fact equal to - T(t), 
we could write its average as 

-E{T(t)\ = -[ Fr{W(u) = 0) du. 
•'o 

But a may be smaller in magnitude than T(t); this fact explains the 
presence of the conditional probability in the integrand of 

E{ a } = - [ Pr{K(t) - K(u) 

- t + u ^ w\ W(a) = 0} Pt{W(u) = 0} du. 

The kernel in this expression is a probability, so it reduces the mag- 
nitude of the integrand whenever it is less than one. 

X. DELAY EXAMPLE: POISSON ARRIVALS, GENERAL SERVICE TIMES 

For a first example, we assume that customers arrive in a Poisson 
process of intensity X, and that service times are mutually independent, 
with a general distribution function B(x). Such a system has been 
treated before, 9 ' 10,11 but few explicit formulas are known except for the 
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Fig. 10 — Graph of a(w). 
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exponential case 

x ^ 0, 

x < 0, 

considered by Ledermann and Reuter 11 and Bailey. 9 The author 10 has 
sketched a method for calculating Pr{W(t) = 0| W(0)} for general 
distributions B(x) of service time, but gave no explicit results. He 
proved that the Laplace transform of Pr{ W(t) = | T'T(O)} is given by 

-,(t)IF(0) 

Pr [W(t) = | 1F(0)| dt = , . , Re(r) > 0, (36) 



/; 



where 77(7-) is the unique root in the right half -plane of the equation 

r - V + X = A£*(7/), Re(r) > 0, 
with 



B*( s ) = r e~ st dB(t). 
Jo 



It was also shown that any function of r?(r), analytic in the right half- 
plane, could be expanded in a Lagrange series. We shall now derive 
these results (by a quite different way) directly from the general in- 
tegral equation (9), and then obtain some specific new results. 

Since arrivals are Poisson and service times are independent with dis- 
tribution B(x), then the load process K(t) is the compound Poisson 
process, and 

jp, -slK(t)-K(u)]i _ \(t-u)[B*U)-l] 

The kernel R(t,u,0) of (9) is 

E " X " ( / ~ ^ B n (t - u) = Pv{K(t) - K(u) ^ t - u], (37) 

«=o nl 

where B„(x) is the convolution of B(x) with itself n times, i.e., the dis- 
tribution of the sum of n service times, and B (x) is 1 for x 22 and 
is otherwise. In fact, the Poisson term in (37) is just the chance that 
n customers arrive in (u,t), and B„(t — u) is the chance that their com- 
bined service time is not more than (t — u). Similarly, we find that 
t 
Py{K{1) ^ u) du = E{max[0,t - K(t)]\ 

= £ e l ^ t)n [ max[0,f - x- W(0)] dB n {x). (38) 

n=n nl Ja 



/. 
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Equation (9) for this example is therefore 

E e ~*'M" f max(0,« - x - TF(0)) dB n (x) = 
"=° n! ■'■ (39) 

I E ^-V — 5.0 - tt)Pr{T7(tt) = | 17(0)} du. 

Since the right-hand side is a convolution, we take Laplace transforms. 
That of the kernel, (37), can be written as 

f C'Rit) dt = ±-.r I* 1 " e - tl *-° +x -' B ' ] * <ft, c > 

Jo 2irl Jo Jc-iv s 

ds (4°) 

i fc+ioo 

= JL( 5 

2iri Jc-ioo t - s + X - \B* ' 

since the order of integration can be interchanged. 

Let S R be the semicircle that is the right-hand half of the circle 
| s — c | = R. It can be seen that on this semicircle 

T _ 5+ 5 ~x- Xfi * = 0(fi - 2) > 

so that 

s~ ds 



«-.«, 2wi Js R T — s + X — XJ3* 






It has been shown 10 that, for Re(r) > 0, the function r — s 4- X 
XB*(s) has a unique zero, 77(7), in the right half-plane. Hence, 



•+ioo -1 



2ri Je— too r — 



S rfs 



2?rt Jc-i=o r - s 4- X - \B* 

(41) 

— ^-r- [residue of (r — s + X — XB*) _1 at s = t/(t)], 
17(7) 

when c < Re tj(t). 

The Laplace transform of E{max[0,t - K(t)]} can be written in the 
form 

— / / <r' (T ~ 8 +X) E — (B*)" ^-=- ds *• ( 42 ) 

2?ri Jo Jc-i» »=o n! s 2 

Formula (42) simplifies to 

2W ic-i- r - s + X - \B* ' 
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which, by arguments like those already used for (40), can be shown to 
equal 

h(r)]"V" (T)w,(0) [residue of (t - s 4- X - X£*) _1 at s = tj(t)]. (43) 

It follows from (39), (40) and (43) that the Laplace transform of 
Pr{T^(0 = | W(0)} is given by 

(■oo -v(t)W(0) 

/ e~ Tt Pr{W(t) = I W(0) } dt = -^— , (44) 

where jj(t) satisfies r — v\ + X = Xf?*(r?). 

As shown, 10 any function F of t?(t), analytic in Re(r) > 0, may be 
expanded in the Lagrange series 

FHr)] = F(r + X) + £ ^ fl \f ( B *r] . 

n =i n\ as" 1 ]_ds J 8 =t+x 

We can use this expansion to invert the transform given by (44). In 
some cases, this inversion gives an explicit expression for 

Vr{W(t) = 0|T^(0)} 

in terms of the kernel and forcing function of the integral equation (39). 
Setting F(x) = af'e - *"'* * in the expansion and inverting the resulting 
transform, we find 

Pr{W(t) = | T^(0)| = e~ u U[t - W(0)) 

g e u x n t n 1 1 r Bn[u _ w(0)] du + w(0)Bn[t _ wm \ m 

»=i n! [Jo J 

By rearranging terms in (45), comparing with (37) and (38), and re- 
calling that TF(O) = K(0) by convention, we can put (45) into the 
form 

Pr{W(t) = 0|TT(0)} = 

T^jmax^ - K(t)]} + 2^) p r |^ (<) ^ ,j. (46) 

Note that E{max[0,t — K(t)]\ is the left-hand side (the forcing func- 
tion) of the integral equation (39) and that, when W(0) = 0, (46) gives 
an explicit representation of the solution of (39) in terms of the forcing 
function alone. The intuitive meaning of (46) can be expressed as fol- 
lows: the chance that the system is empty, conditional on the initial 
load W(0) is equal to 
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average of greater of and 1 — C K(t) 

+ — — (chance that TF(0) plus load arriving in (0,0 is at most /). 

l> 

It is easy to obtain new specific formulas from (46). For example, 
suppose that service times have the fixed length b. In this case of "con- 
stant" service times, for t > TF(0) 

E[max[0,t - K(t)]\ = £ e <( *°" [t - rib - IF(0)J, 

nb^l — JF(0) TCI 

Pr\K(t)£t\- E ^^-", 

n b^l-W{0) n\ 

and hence 

Pr \W(t) = | ir(0)} = 1 - X6 - P(T,\t) + X6P(T - U0, 
where bT = t - TF(0) and 

P(c,a) = E 



is the cumulative term (the "tail") of the Poisson distribution with 
mean a. This formula for Pr{TF(0 = 0| TF(0)} for constant service 
times was used to compute the curve of Fig. 1. 

By rewriting (46) in terms of inversion integrals, we obtain another 
representation of Pr{TF(0 = 0| TF(0)}. This one is more useful be- 
cause, from it, we can find explicit, formulas for new cases (by evaluating 
the inversion integrals). From (42) we see that 

ds 



-«1F(0)+X(B*(«) _ 

s 2 



E {mux [0,1 - K(t)]} = %-. \ ">-" 

Ztn •> c—i oo 

and (45) yields 

pr{*«) % t) = £! pv-^ ^*^-. 

Therefore, from (46), 

Pr{IF(/) =0|1F(0)| ^f'^^ TOTl ^^ (4?) 

As we shall see, this result is useful because it is often easier to evaluate 
the complex integral than to sum the series (37) or (38). 
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For example, if service times have the negative exponential distribu- 
tion with mean u, then it has essentially been shown 9,11 that 

^Vv{W(t) = 1 W(0) = 0} = - (*Y 6 -o.+x>*/i[2(/*)*«] f 



dt l ' W t 

Hence 

Vv{W{t) = 1 17(0) -Oj-l-pf p-*e- ll+p) * h(2p>x) -, (48) 

Jo x 

where p = X//x = traffic intensity, and tit is the time measured in mean 
holding times. This result can be obtained directly from (44) by solving 
a quadratic equation, and it can be put into another form by using (47). 
Thus, 



— (X +/i) t /• c+ii— i oo 

Pr{TF(0 = 0|TF(0) = 0} = ^-— e ul+h " lu (u - M ) 



Writing the integrand as e ut (u — /x) 2 (e x "' /u — 1 + 1) and using Ref. 
12, p. 244, no. 31, we find 

Pr{W(t) = 1 W(0) - 0} = 

[ r t/(4X M ) / ,.2 \ "I 

1 + /„ .-"^ (!-,&)«»)*]■ 



(' 



For another example, suppose that service times have the "gamma" 
probability density 

r(i) 

whose Laplace transform is [n/(n + s)] . Then (47), with a change of 
variable, gives 

— (K+ii)l /•c+ieo 

Pv{W(t) = 0|1F(0) = 0} = - ir — / e «*»'W*« (u - nfdu, 



and 



Pr{W(t) = 0\W(0) = 0} = fV* 1 f e"""G(w)(l + M u) du} , 



t This formula is a simplification of that of Ref. 9 obtained by using standard 
Bessel function relations; /i(x) is the Bessel function of imaginary argument, of 
order one. 
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where 

f,/ \ -5/2 -3/4 / 3/2-i'/bu 2 t/ \ j^ 

G(u) = u a I x e h{x) dx 

Jo 

and a = 4X/i . 

XI. PRECISE DERIVATION OF LOSS RESULTS 

The proof of (10) given in Section VIII is rigorous up to and includ- 
ing (17), so it suffices to give a precise construction of the kernels 
R(t,u), G(t,u) and Q(t,u). Let x n be the nth epoch at which A(t) be- 
comes equal to 0, and y n be the arrival time of the nth successful unit. 
It is readily seen that 

E{6(t,t + w)} =2 Pr {?/« ^ t and Z'H-i ^ * + w \ 

n=l 

DO ft 

= 2 / Pr{z»+i £ t+ w\y n = u] d Pr{?y„ ^ n} . 

n=l Jo 

Similarly, 

E[S(t)} = f)Pr{2/„ ^ *}. (49) 

Thus, each Pr{?/„ ^ • } measure, n = 1, 2, • • • , is absolutely continuous 
with respect to (49). We can therefore express 2?{0(M + w)} as the in- 
tegral of a kernel R(t + w,w) against #{£( •)} measure, as in 

E{6(t,t + w)} = f flU + w,M)dtf{£(n)}, 

Jo 

where the kernel is defined in terms of the indicated Radon-Nikodym 
derivatives by 

B(ftti) = E Pr{x,H4 ^ y|y. = n} ^ {<S(m)} ■ 

In a hke manner, we can write 

Pr{A(0 = and some successful arrival occurred in (0,t]} = 

(50) 

n=l 

Introducing the kernel 

G(y,u) = ± ?v{x n+ i ^ y < y n+l \y n = u] d % { ^" ] , 
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we can render (50) as 

f G(t,u)dE{S(u)}. 
Jo 

The kernel Q(t,u), finally, can be defined as R(t,u) — G(t,u) or as 

Ed r ^ t\ i d Pr{?/„ ^ u) 

The sense in which Q( • ,u) is a "distribution function" is given by the 
following result : for almost all u with respect to E{ S( • ) j measure, Q(- ,u) 
is a distribution function. We show first that the derivatives 



/ % _ dPr{y n ^ y) 
*" (W) ~ dE{S(u)\ 



have the properties 



£*(»)-!. (51) 

n=l 

=g *.(«) ^ 1, (52) 

almost everywhere in E{S( • )}. Now (51) is true by definition. Suppose 
(52) failed on a set B of positive measure; then either E{S( • )} is not a 
positive measure, or else 



7>r[yn*B} > f dE{S(u)}, 

J n 



both of which are impossible. It is readily seen, by an elementary de- 
composition of the events, that, for each n ^ 1 , 

Pr{.T n+ i ^ 1 1 ijn = u] — Pr{.T n+ i ^ t < y„+i\ y n = u) = 

Pr{z/„ + i ^ t\y n = u}. 

Except on a set C n of Pr[y n ^ • } measure zero, this is a distribution 
function in t, and the result follows from 



/. 



dE{S(u)} = 0, if C = U Cn. 



C 71=1 



XII. A LIMIT THEOREM 

We shall prove that dependence of the kernels R(t,u) and Q(t,u) on 
(t — u) only and existence of the "mean service time of successful ar- 
rivals" are sufficient to guarantee that 

Pr{A(0 ^ w] 
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approaches a limit as t — > °° . It is natural to study cases in which only 
difference, kernels occur, because of the Volterra equation derived in 
Section IX. The result to be proved requires no restriction on the "ar- 
rival rate" of units — an upper bound is unnecessary because of the loss 
operation; i.e., if the arrival rate increases, the rate at which successful 
units leave the system can only increase to a limit. 

Theorem: If the kernels R(t,u) and Q(t,u) only depend on (t — u), 
if the average 

6 = / [1 — R(x)]dx = (average service time of a successful unit) 
Jo 

exists, and if Q( ■ ) is not a lattice distribution, then 

a 4- [1 - R(u)] du 
Jo 



where 



limPr{A(0 ^ w) = 

t-"» a -r o 



(a+b) = [ [1 - Q(x)]dx, 

Jn 



(53) 



the limit being one if (a + b) = ». 

Proof: By the remarks at the end of Section XI, both R(x) and Q(x) 
may be taken to be the distribution functions of positive variates. Equa- 
tion (11) becomes a renewal equation, and E{S(u)} is essentially the 
renewal function H{ • ) of Smith. 13 The integrand of ( 10) is 

1 - R(t + w - u), 

and is nonincreasing and integrable. Also, Q( ■ ) is not of lattice type. 
Hence, by Smith's Theorem 1 : 

pt [ [1 - R(u)] du 
lim / [1 - R(t +io- u)] dElS(u)} = -^—- — , (54) 

/ xdQ(x) 

Jo 

which gives (53) upon rearrangement if (54) is taken to be zero if the 
mean of Q(x) does not exist. 
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